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ABSTRACT 

Hubble Space Telescope images of the galaxy cluster Abell 2261, obtained as part of the Cluster 
Lansing And Supernova survey with Hubble, show that the brightest galaxy in the cluster, A2261- 
BCG, has the largest core yet detected in any galaxy. The cusp radius of A2261-BCG is 3.2 kpc, twice 
as big as the next largest core known, and ^ 3x bigger than those typically seen in the most luminous 
BCGs. The morphology of the core in A2261-BCG is also unusual, having a completely flat interior 
surface brightness profile, rather than the typical shallow cusp rising into the center. This implies that 
the galaxy has a core with constant or even centrally decreasing stellar density. Interpretation of the 
core as an end product of the "scouring" action of a binary supermassive black hole implies a total 
black hole mass ~ IQ^^Mq from the extrapolation of most relationships between core structure and 
black hole mass. The core falls Icr above the cusp-radius versus galaxy luminosity relation. Its large 
size in real terms, and the extremely large black hole mass required to generate it, raise the possibility 
that the core has been enlarged by additional processes, such as the ejection of the black holes that 
originally generated the core. The flat central stellar density profile is consistent with this hypothesis. 
The core is also displaced by 0.7 kpc from the center of the surrounding envelope, consistent with a 
local dynamical perturbation of the core. 

Subject headings: galaxies: nuclei — galaxies: photometry — galaxies: structure 



1. A LARGE CORE AS A TEST OF CORE FORMATION 

Brightest cluster galaxies and most luminous early- 
type galaxies brighter than My ^ —21 have "cores" in 



their central starlight distributions (Faber et al. 1997 



Laine et al. 112002 Lauer et al. 2007a I. Cores are markec 



by a distinct physical radius interior to which the pro- 
jected starlight surface-brightness increases only slowly 
as r — ^ 0, in marked contrast to the surrounding enve- 
lope, which has a much steeper profile (in logarithmic 
units). In qualitative terms, a core looks like a central 
"plateau" in the starlight distribution. In more quanti- 
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tative terms, a core can be defined as the central region 
of a galaxy where the surface brightness takes the form 



of a shallow cusp, / 



(Lauer et al 
than 



"MI 



1995 



2005) 



cx , with 7 < 0.3 as r — >■ 
Importantly, galaxies fainter 
having 



generally do not have cores 
7 > 0.5 instead, as r — 0. This distinction is of physical 
interest, as the presence or absence of a core correlates 
with the strength of the stellar rotation field, isophote 
shape, nuclear radio emission, and overall X-r ay emis 



sion, in addition to the total galaxy luminosity (Faber et 
al. 1997| [Lauer eral]|2007b[ ) 



The lormation of cores has long been thought to be 
due to the action of black holes on the central struc- 
ture of galaxies. Their form and size may refiect both 
the mass of the central black hole in the gal axies, and 
the rnerger h istory that created the galaxies. |Begelmaii| 



et al. ( 1980 1 hypothesized that a binary black hole cre- 
ated in the merger of two galaxies would eject stars from 
the center of the newly created system as the binary 
slowly hardened. In simple terms, the black hole binary 
"scours" out the center of the galaxy, thus "flattening" 
the otherwise steeply rising central starlight distribution 
as r — )■ 0. Subsequent N-body sim ulations have demon- 
strated this phenomenon, directly (Ebisiiz aki et al.|1991 



stratea tn is pnenomenon , directly (tbisuzaRi 
[Maki no J9971 [Milosavljevic fc Merritt 20dT| T 



Faber et al. ( |1997p showed that cores occur in the most 
luminous elliptical galaxies, and are correlated with slow- 
rotation and "boxy" isophotes in these galaxies, conclud- 
ing that core formation is a natural end-point of dissipa- 
tionless mergers of two progenitor elliptical galaxies. The 
conclusion that ne arly every ellipt ical galaxy has a black 



hole at its center ( Magorrian et al.j|1998j ), coupled with 
the conclusion that the most massive elliptical galaxies 
were formed by merging pre-existing gas-free galaxies, 
explains why cores are found in nearly all luminous el- 
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lipticals. 

While "core scouring" is an attractive hypothesis for 
the formation of cores, there may be additio nal mecha- 
nisms fo r binary black holes to generate cores. I 



&: Reesj ( |1989p suggested that when the two black holes 
in the binary ultimately merge, asymmetric emission of 
gravitational radiation could eject the merged hole from 
the center of the galaxy, causing the center to "rebound" 



ritt et al. (2004), Boylan 


-Kolchin et al. 


( 


2004 


), 


landris &: Merritt 


( 


2008 


) studied this problem 



jua- 



tion of the merged black hole indeed could cause the 
central distribution of starlight to re-adjust such that it 
would create a core in the projected stellar surface bright- 
ness profile. An interesting ancillary effect discussed in 
these works is the possibility that the ejected black hole 
would remain bound to the host galaxy on a radial orbit. 
In that case, it would repeatedly fall through the cen- 
ter of the galaxy, continuing to enlarge the core through 
dynamical friction. 

Some observational support for the scouring origin of 
cores comes from the measurement of core "mass deficits" 
compared to the b lack hole masses in the same galaxies 
( |Faber et al.(l997l ). The mass deficit, Md, is the inferred 



mass ot stars ejected from the center of the galaxy re- 
quited to create a core, and is estimated by a reference 
to a postulated initial form of t he g alaxy fitted to the 
envelope, such as a Sersic law fGraham 20 041). Vari- 
ous studies estimating mass deficits ti-'aber eTal. 111997 



Milosavljevic et al.|2002 


IRavindranath et al.l2U02| |Gra- 


ham 2004 


|Merritt||2006 


Lauer et al.| 2007a|]Kormendy 


& Bender 


2009) typically find Md(x M,, the black hole 



but with large scatter about the mean relation. 

The observational context for understanding the large 
core in the brightest cluster galaxy in Abell 2261 (here- 
after referred to as A2261-BCG) is provided by the exten- 
sive work done with the Hubble Sp ace Telescope (HST 
on the structure of nearby galaxies. Lauer et al.] (j2007al 



constructed a large sample of galaxies that had high- 
resolution surface photometry obtained with HST; sig- 
nificantly, it inclu ded the core parameters measured by 



Laine et al. ( 2002 ), who studied a large sample of bright 
est cluster galaxies. Very few galax ies have cores as large 



as 1 kpc; the largest core in the Lauer et al. (2007a I 
sample is that for NGC 6166 = A2199-BC(i, wh i cF'Eas 
a core size of ~ 1.5 kpc. McNamara et al. (20091 



drew attention to the large core that they discovered m 
BCG MS0735. 6+7421; however, it is compatible w i th the 



largest BGG cores measured in the Laine et al. (20021 
sample. The core in A2261-BCG, however, is over twice 
as large as that in NGC 6166. It provides an extreme 
test of the possible mechanisms hypothesized to gener- 
ate cores. 

2. OBSERVATIONS AND PROPERTIES OF A2261-BCG 

The center of Abell 2261 was observed for a total of 20 
orbits as part of the Cluster Lensing And Supernova sur- 
vey with Hubble (CLASH) multi-cycle treasury program 
between March 9, 2011 and May 21, 2011 in 16 jjroad - 
band filters from 0.22 — 1.6 /xm ( Postman et al.|[2012 ). 
For the calculation of physical quantities, we assume a 
cosmology with ^lm — 0.3, f^A = 0.7, and Hq = 70 



km s"^ Mp c~^ At t he mean redshift of Abell 2261, 
z = 0.2248 ( |Coe et al.|[2012j ), 1" subtends 3.61 kpc and 
the distance modulus is 40.241. 

AbeU 2261 is included in the CLASH X-ray selected 
subsample of 20 clusters. The CLASH X-ray selected 
sample consists of clusters with X-ray temperatures 
greater than 5 keV and exhibit a high degree of dynam- 
ical relaxation as evidenced by Chandra X-ray Observa- 
tory images that show well-defined central surface bright- 
ness peaks and nearly concentric isophotes. The intr- 
acluster medium (ICM) of Abell 2261, in particular, is 
characterized by an X-ray temperature of T^ = 7.6±0.30 
keV, a bolometric X-ray luminosity of 1.80 ± 0.20 x 10^^ 
erg s~^ and a [Fe/H] abundance ratio that is 0.31 ± 0.06 
times the solar value ( Anders fc Grevesse|1989 ). The esti- 
mated mass within r25oo (radius where the density is 2500 
times the critical density) is M2500 = (2.9±0.5) x 10^'* M© 
w ith a gas fraction of .115 ± 0.01. 



Maughan et al.| ( |2008[ ) found Abell 2261 to have a small 
level ot substru cture m its X -ray gas surface brightness 



distribution and Gilmour et a l. ( 2009 ) classified the clus 
ter as "disturbecl?^ Coe et al. (2012) find that applica tion 
of the caustic tech nique (e.g., Diafcrio fc Geller,1997 Di- 
aferio et al.[2005 ) to spectroscopically measured galaxies 
in the vicinity ot Abell 2261 suggests that the dynamical 
center of the cluster is located ^ 6' (~1.3 Mpc) south of 
the BCG. However, the BCG (z = 0.2233) in Abell 2261 
is at equatorial coordinates of 17:22:27.18 +32:07:57.1 
(J2000), putting it within i;'6 (5.8 kpc) of the centroid of 
the ICM X-ray emission. Its mean velocity offset relative 
to that of the cluster mean redshift is 367 km s~^. The 
distribution of BCG velocity offsets for a sa mple of 42 
Abell clusters has a mean value of 264 km s^^ (Postman 
An expanded dataset of 174 Abell clus 



& Lauer 



1995) 



ters, each with at least 50 spectroscopically confirmed 
members, yields a mean value of the BCG velocity offset 
of 234 km s~^, with 22% having offsets of at least 350 
km s^^. The A2261-BCG, thus, appears to be reasonably 
aligned with the center of the cluster's main gravitational 
potential well. 

The unusual core of A2261-BCG was discovered dur- 
ing the initial inspection of the CLASH HST images of 
A2261. Visually, the core presents itself as a large, round, 
uniform disk of low surface brightness, as can be seen 
in Figure IT] which shows a color composite of the cen- 
tral 2' X 2^region of AbcU 2261 made from the CLASH 
images. In addition, four compact sources with colors 
similar to the BCG itself are superimposed on the out- 
skirts of the core. The 3 brightest knots are marginally 
detected (< 3cr) in the WFC3/UVIS F336W image, well 
detected (> 6cr) at longer wavelengths, and not detected 
in either the F225W or F275W UVIS passbands. 

2.1. Nuclear Activity 

Understanding the central structure of A2261-BCG re- 
quires knowledge of whether or not the core hosts a cen- 
tral supermassive black hole. This problem will be con- 
sidered at length later in the paper, but we note here that 
the center of the galaxy does harbor a radio source that is 
at least an order of magnitude more powerful than those 
associated with star forming regions, and is thus evidence 
in favor of an active galacti c nucleus (AGN). Th e NRAO 



VLA Sky Survey (NVSS; 'Condon et aL 



source at 1.4 GHz with an integrated Mux o: 



1998) finds a 
^:^0.5 mJy 
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Figure 1. Color composite HST image, from CLASH ACS/WFC and WFC3/IR images, showing the BCG in A2261 and its neighbors 
in the central 2x2 arcminute region of the cluster. The insert in the lower right hand corner shows a zoomed in region centered on the 
BCG with contrast adjusted to highlight the bright knots (labelled 1,2,3,4) in the core. The orientation is north up and west to the right. 
The faint "figure 8" patterns at the 6 o'clock and 11 o'clock positions are due to internal reflections in the ACS camera of light from a 
nearby bright star. The red "diamond" at the 10 o'clock position near the BCG is caused by a gap in areal coverage due to the multiple 
orientations used in the CLASH survey. The red "blob" at the right edge of the image is a WFC3/IR detector artifact that does not easily 
calibrate out. 

resolution is 45"). 

If the radio source is at the cluster redshift then the 
FIRST detection yields an absolute luminosity of 7.8 x 
Iq39 gj.g g-i 1.4GHz corresponding to I^^Hz = 5 x 
1023 w Hz-i. [Hlavacek-Larrondo et al] ( [20121 also find 
a central source with luminosity 5.0 x 10^^ erg s~^ at 
5 GHz corresponding to L^ qghz = 1 x lO^^ W Hz"^. 
The radio fluxes are indicative of a mildly powe rful FR I 
ra dio sou rce typical of those seen in BCGs (e.g. Croft et 



that lies 6 "3 east of the BCG position. No other NVSS 
sources are found within 2 ^ 45 of the BCG. Data from the 
Faint Images of the Radio Sky at T wenty-Centimeters 
(FIRST) survey ( |White et aL||l997[ ) reveals a 3.4 mJy 
source that is more comcident with the core, lying 1''6 
west of the BCG position. The FIRST detection limit 
at the source position is 0.99 mJy per beam. No other 
FIRST sources are found within l! 92 of the BCG. Pre- 
sumably the NVSS and FIRST detections closest to the 
BCG correspond to the same source. The two radio posi- 
tions are 7 "5 apart, which is less than a 2a difference and 
is consistent with the different astrometric uncertainties 
for the two surveys. Unfortunately, the VLA sky survey 
data currently do not have sufficient angular resolution 
to determine whether one of the knots, or the core itself, 
may host the AGN (VLA FIRST resolution is 5", NVSS 



al. 



2007) 



Apart from the radio source, there is also a faint 24 jim 
detection of 0.44 mJy (Hoffer et al. 20121, which would 
indicate either star forniation or a hot dust torus around 
an AGN. As in the case of the radio observations, how- 
ever, the angular resolution is probably too poor to dis- 
tinguish whether or not the flux comes from the center 
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Figure 2. Spectrum of the BCG in Abell 2261. Its redshift is 
z = 0.2232. This spectrum is from the GMOS spectrograph on 
the Gemini-N 8-meter telescope (P.I. H. Hockstra). Prominent 
absorption features are noted. The data from the GMOS blue and 
red channels are denoted, respectively, by the blue and red traces. 
The break at 5100A is due to a small data gap between the two 
channels. This spectrum was used to derive the stellar velocity 
dispersion of the BCG reported in this work, a = 387 it 16 km . 

of the core or any of the associated knots. 

In contrast, optical spectra do not reveal any signifi- 
cant evidence of emission lines that may be associated 
with an AGN or cooling flow activity (the spectra are 
presented in the next subsection). Further, there are 
no known X-ray or EUV sources identified with any of 
the compact knots; the X-ray limit on any p oint source 



in the core is I/2-ioke V < 3.8 x 10*^ erg s ^ (|Hlavacek- 



Larrondo et al. 2012). At the same time, it is typical 



for central radio source in BCGs not to have strong X- 
ray counterparts. A black hole mass of ~ 10^° Mq is in 
general consistent with the power requirements associ- 
ated with the ICM shock fronts seen in other clusters 



(e.g. McNamara et al. 2009); but for such high mass 
black holes, the power needed for such shocks is only 
^ 0.001 L/Led, which therefore im plies that the black 
holes can b e radiatively inefficient. |Hlavacek-Larrondo, 
et al. (2012) note that synchrotron cooling can strongly 
attect the X-ray luminosity from the most massive black 
holes since the cooling break occurs below X-ray wave- 
lengths, making such massive black holes under luminous 
at X-ray wavelengths compare d to their radio luminosity. 



For A 2261-BCG in particular, Hlavacek-Larrondo et al. 
(2012) estimate an Eddingtion ratio of 



10- 



10- 



2.2. Photometry 
The details of the image reduction and co-addition of 



et al. 


(2012 


) and 


Coe et al. 


(2012 



HST CLASH images. The ACS photometry is on the 
AB magnitude system but we will express much of the 
reduced photometry in r est-frame V-band (V ega-based) 



for compatibility with the Lauer et al. ( 2007a ) analysis of 
central galaxy photometry. We use a suite of Bruzual ^ 
Chariot (2003 hereafter BC03) synthetic stellar popula- 
tion models to derive linear photometric transformations 



from the observed, extinction-corrected ACS / WFC mag- 
nitudes and colors to the rest-frame Johnson V-band, Vo ■ 
The BC03 models used to empirically derive the trans- 
formation equations have the following parameters: ex- 
ponential star formation rate e-folding times, r, of Gyr 
(SSP model), 0.2 Gyr, 0.6 Gyr, and 1.0 Gyr; metalhci- 
ties of 0.25^0, I.OZq, and 2.5^0, and ages from 2 Gyr 
to the age of the universe at the cluster redshift in 0.5 
Gyr intervals. These models more than span the range 
of observed cluster galaxy SEDs. To compute photo- 
metric transformations, we take the BC03 spectral en- 
ergy distributions (SEDs) for each of the above models 
and compute the observed photometry and colors at the 
cluster redshift and then use the same SEDs to compute 
the rest-frame Johnson V magnitude. The transforma- 
tion equation parameters are then derived using a linear 
least-squares fitting procedure. In addition, we compute 
an evolution correction, Ce„, from a BC03 t — 0.6 Gyr so- 
lar metallicity model with a formation epoch of z = 4.5. 
The transformation equation is: 



Vo = (F814W + Ce,,F814W 

F814W 



0.5186 X (F606W- 

Cev,F606W — Cev.FSUw) + 0.1764 

(F814W) -1-0.5186 X (F606W-F814W)+ 0.4213 (1) 



where F606W and F814W are the Galactic extinction- 
corrected ACS/WFC measurements in AB-magnitudes, 
Cev,F606W = 0.256, Cev,F8iiW = 0.233, and Vo is on 
the Vega system. The extinction corrections used are 
0.127 mag and 0.079 mag for the F606W and F814W 
filters, respectively. The rms scatter about the transfor- 
mation equation shown in equation 1 is 0.006 mag. We 
tested the robustness of the method by computing simi- 
lar transformations for F625W and (F625W-F814W) and 
for F850LP and (F814W-F850LP). They yield the same 
rest-frame Vo values to within ±0.01 magnitudes. 

The total luminosity of A2261-BCG was determined by 
fitting its surface photometry with an r^^^ la w. This is 



consistent with the methodology in Lauer et al. (2007a 
which will be used to provide the context for the present 
galaxy. The BCG is highly luminous with a total abso- 
lute magnitude of —24.70 in the V-band (see Table m). 
Because we want to compare this to the z = sample 
of BCGs, this luminosity includes an evolutionary cor- 
rection of -1-0.26 mag to account for the aging of the 
stellar population since the z — 0.22 epoch. Even then, 
A2261-BCG is among th e most luminous BCGs known 
( [Postman k Lauer||1995 ). 

bom sources tall within the outskirts of the core. Their 
coordinates and photometric properties are summarized 
in Table [T] The brightest source, "knot 3," is well re- 
solved. Its profile was m easu red simultaneously with the 
BCG (as is described in { 3.1 ) and has a roughly exponen- 
tial form. The close pair or compact sources north-east 
of the core center, knots 1 and 2, are both marginally 
resolved with 0''05 half- intensity radii. The faint source 
south of the core, knot 4, is unresolved. The photometric 
errors in Table [T] given here include both the random and 
estimated systematic errors. The random errors (photon 
shot noise, read noise, dark current) are in the range 
0.005 mag to 0.010 mag. The systematic errors are rela- 
tively small but still about 2 to 4 times larger than the 
corresponding random errors. This is due to spatial vari- 



The Core of A2261-BCG 



5 



ations in the background level after the BCG model is 
subtracted from the image. To estimate the amplitude 
of the systematic error we compared the knot photom- 
etry derived from the BCG-subtracted image with that 
derived by using a local sky subtraction (instead of BCG 
subtraction) where the background is estimated from an 
annulus with an inner radius of 0"6 and an outer radius 
of 1''2 centered on each knot. The two measures of the 
photometry typically agree to within 0.015 mag to 0.020 
mag. Based on this, we adopt 0.018 mag as an estimate 
of the systematic error. 

2.3. Spectroscopy and The Central BCG Velocity 
Dispersion 

Moderate resolution spectroscopy of the BCG exists 
from the Sloan Digital Sky Survey (SDSS) DR8 and from 
the Gemini Observatory data archive. The Gemini data 
are deep long-slit spectra with GMOS-N from the pro- 
gram GN-2008A-Q-103 (P.I. Chris Bildfell). The AbeU 
2261 observations were obtained on March 16, 2008. The 
GMOS-N long-slit, with 0"75 width, was placed on the 
BCG and oriented at a position angle that intersects two 
of the three bright knots in the northern part of the core 
(but not the brightest one). The spectra are split into 
dual red and blue channels to optimize sensitivity. The 
full spectral range is ~ 4140 — 6040A with a spectral res- 
olution of i? ~' 1000 at the blue end. The ID co-added 
CMOS spectrum is shown in Figure [2] 

We measured the stellar velocity dispersion (cr) within 
the core of A2261-BCG from the CMOS spectrum using 
the ID L code pPXF developed by |Cappellari fc Emsellern 
(2004), which fits a linear combination of template spec- 



tra to the observed galaxy spectrum to minimize tem- 
pla te mismatch. The tem plate spectra are provided by 
the Vazdekis et al. (20101 single burst models for a range 
of ages and metallicities, which are based on the empiri - 
cal MILES stellar library fSanchez-B lazquez et al.|2006 ) . 

The fits are computed separately in the red and blue 
channel spectra, covering rest-frame wavelengths 3700- 
4150A and 4250-4700A, respectively in order to avoid 
regions near the chip edges where the wavelength solu- 
tion is less reliable. The bright skyline at the 5682-5688A 
Na doublet is masked in the fitting procedure, and the 
fits are computed weighting all pixels equally. We find 
consistent values of cr = 393 ± 13 km s~^ on the blue side 
and (7 = 380±8 km s~^ on the red side. The higher-order 
non-Gaussian velocity moments H3 and HA derived from 
the fits are negligible. We therefore adopt a = 387 ± 16 
km s^^ as an estimate of the ensemble stellar velocity 
dispersion, including systematic errors. The Gemini ve- 
locity dispersion value is in excellent agreement with the 
388 ± 19 km s~^ velocity dispersion estimate from the 
DR8 SDSS database. This is among the highest central 
velocity dispersion values known. 

3. THE CENTRAL STRUCTURE OF A2261-BCG 

3.1. The Morphology of the Core 

The surface brightness profile of A2261-BCG was mea- 
su red from the F814W imag:e deco nvolved with 20 cycles 
of Lucy ( 1974| ) - Richardson] ( 1972|) deconvolution to cor- 
rect tor the blurring of the FSl*'. I'he deconvolved image 
of the core is given in Figure [Sj Deconvolution works well 
on HST images for recovering estimates of the intrinsic 



light distributions of galaxies (Lauer et al. 1998, 2005). 
Given the angularly-large and hat-protile or the A226: - 
BCG core, the effects of deconvolution in the present con- 
text are extremely modest. The increase in the central 
surface brightness within the core after deconvolution is 
only ~ 3%, with similar effect on the measured angu- 
lar size of the core. The real import of deconvolution 
for A2261-BCG is to recover the forms of the compact 
sources falling around the core, and to reduce the effects 
of their scattering "wings" on the core, itself. 

Figure [3] also shows the residuals obtained after a 2D 
model reconstructed from the brightness profile of the 
core was subtracted from the ima ge. The prof ile for 
r > 0"5 was measured using the Lauer (1986) algo- 
rithm (operating in the xvista image processing sys- 
tem), which solves for the overlapping light distribu- 
ti ons of multipl e galaxies; the high-resolution algorithm 
of Lauer ( 1985 ) was used interior to this. The surface 
brightness profile is presented in Figure [5j The profile 
ratifies the visual impression that the core is essentially 
flat. In contrast to the typical centr al structure of mos t 
elliptical galaxies (but not all — see Lauer et al.||2002" ), 
there is no sign of any rising cusp in starlight as r — > 0. 
The lack of a cusp actually made it difficult to determine 
the precise center of the galaxy. We did this by taking 
an intensity centroid over the entire core. 

At radii well outside the core, HST-ha,sed surface 
brightness profiles become vulnerable to systematic sky 
measurement errors and so on, due to the low intensity 
levels of the outer isophotes, the angularly-small pixels, 
and limited fields of the cameras. To better character- 
ize the envelope of A2261-BCG, we augmented the ACS 
photometry with surface photometry measured from a 
ground-based R-band image obtained with the Suprime- 
cam imager on the Subaru 8-m telescope. The Subaru- 
derived profile was scaled to match the ACS profile over 
5" < r < 8", then blended with it over the same interval, 
and used solely as the brightness distribution for r > 8". 
The Subaru profile extends to ~ 24" and equivalent V- 
band surface brightnesses fainter than 25 mag/arcsec^. 

Despite the smooth structure of the profile, the resid- 
uals within the core show an number of interesting fea- 
tures. In addition to the presence of the four sources 
noted earlier, the residuals show a dipole-like pattern 
outside the central flat portion of the profile. The sim- 
plest interpretation is that the core is slightly displaced 
from the center of the surrounding envelope by ^ 0"2 
(0.7 kpc) towards position angle ~ 300°. This conclusion 
is supported by the contour map of the core, shown in 
Figure [4j The contour lines outside of the core to the 
SW are clearly spaced closer together than those in the 
NE direction. 

3.2. Analysis of the Surface Brightness Profile 



We de scribe the profile with a "Nuker-law" (Lauer et 
al] [l995| . 



/(r) = 2(^-^)/"/6(^)' 









y 


1 + 











(7-/3)/a 



(2) 



which models the surface-brightness distribution as a 
"broken" power-law. The envelope profile has the form 
I{r) ^ r~^ as r — oo, while the inner cusp has I{r) ^ 
r~^ as r — > 0, with the transition radial-scale provided 
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Table 1 

Abell 2261 BCG and Bright Knot Astrometry and Photometry 



Object 


R.A. 


Dec. 


F814W 


F606W-F814W 


V Mag 


Mv 


ID 








TAR mao-'l'^ 


V vega mag; 


I vega mag; 


BCG 


17:22:27.18 


+32:07:57.30 


14.69 


0.82 


15.53 


-24.70 








(< 0.01) 


(< 0.01) 


(0.01) 


(0.02) 


Knot 1 


17:22:27.23 


+32:07:57.65 


21.53 


0.86 


22.15 


-17.85 








(0.02) 


(0.04) 


(0.06) 


(0.06) 


Knot 2 


17:22:27.21 


+32:07:57.56 


21.89 


0.85 


22.50 


-17.49 








(0.02) 


(0.04) 


(0.06) 


(0.06) 


Knot 3 


17:22:27.13 


+32:07:57.59 


20.20 


0.83 


20.80 


-19.20 








(0.02) 


(0.04) 


(0.06) 


(0.06) 


Knot 4 


17:22:27.14 


+32:07:55.85 


22.93 


0.97 


23.61 


-16.63 








(0.02) 


(0.04) 


(0.06) 


(0.06) 



^ Corrected only for Galactic Extinction. 

^ Corrected for Galactic extinction and K-dimniing (assuming sources arc at tlic cluster redshift). 

^ Corrected for Galactic extinction, K-dimming, and passive evolution from z — 0.224 to 2: — 0. Sec text for 

details 




sample of core galaxies presented in Lauer et al. ( 2005 1 



Figure 3. The left panel showrs the center of the F814W image of the BCG after |Luc"y| | |1974| . |Richardson| ( |1972| l deconvolution. The 
region shown is 12" X 12"; (43.2 X 43.2 kpc) the intensity scale is logarithmic. Nortn is at tne top and east to tne left. The right panel 
shows the residuals after subtraction of a model reconstructed from the surface photometry of the BCG. The over all structure of the 
residuals is a dipole pattern of positive residuals NE of the core and negative residuals to the SW. This suggests that the core is slightly 
displaced from the surrounding envelope in the SW direction. 

by the "break-radius," r?,. The "speed" of transition be- 
tween the envelope and inner cusp is provided by a, while 
lb, the surface-brightness at rt, gives the overall intensity 
normalization. 

The best-fitting Nuker-law is plotted in Figure [5j and 
has the parameters 7 = -0.01, /3 = 1.56, a = 2.41±0.18, 
n = 1.'20, and h of 19.72 in F814W(AB). The fit was 
conducted for r < 8", beyond which the the envelope 
falls away slightly faster than a pure power-law with ra- 
dius. As can be seen, the fit over this domain is excellent, 
with an rms residual of 0.02 mag. The parameters are not 
strongly dependent on the domain of the profile fitted, in 
any case. For a fit conducted limited to r < 3 "7, or half 
of the nominal domain, — 1"15 is recovered, only a 4% 
decrease from the nominal value. The break radius, Vh, 
corresponds to the point of maximum logarithmic curva- 
ture, and thus can also be estim ated by non-parametric 
methods. Lauer et al. (2007a) did this for a the large 



finding the Nuker-fit and non-parametric measurements 
to agree well over a large range in angular core-size, with 
no biases. 

The negative 7 actually implies that the surface bright- 
ness decreases slightly as r — >■ 0. We also performed a 
Nuker-fit forcing 7 = 0, which is essentially identical to 
the nominal fit, but for the central point. In this case, we 
recover ^ = 1.55, a = 2.50, n = 1"30, and h of 19.72. 
The 7 = model falls right on the upper error-bar of 
the central point and thus coincidentally serves as a Icr 
confidence bound on the central profile. 

The center of A2261-BCG may have a luminosity den- 
sity profile that actually decreases as r — >■ 0. This occurs 



in a number of galaxies, including a few BCGs (Lauer et 
ar]|2002[ |2005[). An Abel inversion of the Nuker-modei 
tor A226l-B(JG with 7 = -0.01 indeed shows that r = 
corresponds to a local minimum in stellar density (Figure 
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Semimajor Axis (kpc) 



Figure 4. A contour plot of the core of A2261-BCG. The contour 
levels have an arbitrary zeropoint, but are spaced by 0.25 mag in 
surface brightness. North is to the top and east to the left. Note 
that the contour levels are closer together in the SW direction 
outside the core than they are in the NE, supporting the dipole- 
like residual pattern seen in Figure [3] and the conclusion that the 
core is displaced to the SW relative to the envelope center. 



Semimajor Axis (kpcj 




Semimajor Axis (arcsec) 

Figure 5. The central surface brightness profile of A2261-BCG 
as measured (solid points) is shown with two "Nuker-law" profile 
fits fLauer et al.ll995,). The error bars are smaller than the points, 
but for the central tew measurements. For comparison to previous 
studies the profile is normalized to 2: = V-band (Vega). The solid 
line is the best-fitting Nuker-law and features a slightly depressed 
(7 = —0.01) cusp as r — > 0. The dotted line is an is an r^/'*-law 
(an n = 4 Sersic-law) fitted to the envelope. The triangles indicate 
the cusp-radius. 

[6]). Even inversion of the 7 = model of the brightness 
profile formally implies a ~ 4% decrease in luminosity 
density at the HST resolution limit. This can be under- 
stood as a consequence of the lack of a cusp in combi- 
nation with a sharply-defined core. With a generic sur- 
face brightness profile of the form, I{r) cx (1 -|- 
(both a and /3 positive), the derivative of the profile has 
the form I(r)' cx r" for r « 1. When used with the 
standard Abel-transform, this implies a formal density 




Semimajor Axis (arcsec) 

Figure 6. The implied luminosity density profile of A2261-BCG. 
The density profile resulting from Abel inversion of the best-fitting 
Nuker-law, which features a slightly depressed (7 = —0.01) cusp 
as r — > 0, is shown as a solid line. The density profile plotted as a 
dotted-line is based on inversion of a Nuker-law fitted with 7 = 0. 
The small difference between the two profiles is greatly amplified 
in density as r — > 0. Even the density profile with 7 = is formally 
very slighted depressed as r — 0. The difference between the two 
models is not significant. 



of zero at r = 0, when a > 2, a condition satisfied for 
the A2261-BCG brightness profile. 

Given, however, the modest central decrease in density 
implied by even the 7 — —0.01 model, which just occurs 
near the resolution limit, it is not possible to say with 
confidence if A2261-BCG is yet another example of a 
"hollow core." At the same time there is no sign of any 
centrally rising cusp in central surface brightness, and 
an even constant density core over nearly two decades 
in radius is unusual (see the collection of core density 
profiles in Lauer et"aL]|2007b ). 



3.3. The A2261-BCG Core Compared to Those of 
Other Early-Type Calaxies 

Figure [7] shows that the core of A2261-BCG is the 
largest core yet seen among extensive surveys o f local 
galaxies. The figure is an adaptation of a figure in Lauer 
et al. ( 2007a I , which shows the relationship between the 
core cusp-radius," r^, and total V-band luminosity, Ly 
for a composite sample of HS T studies of the central 
stru ct ure o f ear l y-type galaxi es ([ Lauer et~ al. 1995; Faber] 
et al i997;'Quiricn'ct al. 2000; Ravindranath ct al. 200_3, 
iRest ct al. 2001, Laine et al. 2002; Lauer et al. 2005) 



ihe Laine et al. (2002) study, m particular, tocussed on 
the central structure of BCGs, and is thus particularly 
useful for placing the core of A2261-BCG in context. 

The cusp-radius is the angular or physical scale at 
which the local logarithmic slope of the surface bright- 
ness profile reaches a value of —1/2, as the profile transi- 
tions from the steep envelope to the shallower inner-cusp. 



The V-band used in [Lauer et al. ] l [2007a| and cited in the 
remainder of the paper is Vega- based. 
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In terms of the Nuker-law parameters, 



1/2-7 
/5-1/2 



I/O 



(3) 



This scale was first introduce d by Carollo et al.| (|1997| 
and adopted in the analysis of Lauer et al. ( 2(l07a ) , who 
demonstrated that it gave tighter relationships between 
the core and global properties of the galaxies, than the 
direct use of the break-radius, rj,, did. 
For A2261-BCG, the 7 = -0.01 Nuker-law fitted yields 
= 0"89 ± 0"02, or 3.2 ± 0.1 kpc at the distance of 
A2261-BCG (the 7 = fit also gives = 0"89). As a 
check, this value is consistent with a simple estimate of 
= 0''90, made by fitting a parabola to the core over 
0"7 < r < 1''4. The corresponding "cusp-brightness," 
the surface brightness at r^, is 19.52 mag/arcsec^ 
(F814W) or 20.29 in rest-frame F-band corrected for in- 
terstellar extinction. As shown in Figure[7j this r.y is over 
twice as larg e as the largest BCG cores seen in the |Laine| 
et al. (2002) sample. It is also three times lar ger than the 
large core m MS0735.6-h7421-BCG noted by |McNamara| 
et al.| (|2009|). The A2261- BCG core falls on the high-side 
of the Lauer et al. (2007a) relation between core-size and 
total galaxy-luminosity, 

log(r^/pc) = (1.32 ± 0.11) (-A//v-23)/2.5 + 2.28 ±0.04. 

(4) 

The scatter in log(r^) about this relation is 0.35 dex; 
A2261-BCG falls 0.33 dex above the relation, corre- 
sponding to a l.Ocr deviation. 

4. CONJECTURES ABOUT SUPERMASSIVE BLACK 
HOLES IN THE CORE 

4.1. Estimates of the Black Hole Mass From the Core 

Size 

If cores are created by the scouring action of a binary 
black hole, and the core size is indicative of the mass of 
the binary, then the implied bl ack hole mass for A2261- 
BCG must be extremely large. McConnell et al. (20111 



recently identified two black holes with masses 10 M0 
or greater in the BGGs A1367-BGG=NGC 3842 and 
A1656-BCG=NGC 4889, which ha ve cusp-radii of onl y 
0.3 kpc and 0.7 kpc, respectively (Lauer et al. 2007a), 
implying that M, in A2261-BCG is yet larger than that 
found in those two galaxies. 

The M, — L and Af, — a relations b etween nuclear 
black hole mass and galaxy luminosity (Dressier 1989| 



Kormendy & Richs tone^l995j .M agorrian ct al. 1998) or 



velocity dis persion ( I'terrarese Ife Merritt 2000 ; [{iebhardt 
et al. [20(30 1 give sonie guidance on what to expect tor the 
central black hole mass in A2261-BCG, completely apart 
from its core structure. The M, — L and M, — a (elliptical 
galaxy on ly sample) relations from the comprehensive 



analysis of Giiltekin et al. ( 2009 1 predict M, = 6x10^ 



and 2 x lO'^ Mf?), respectively. The analogous McGonnell 
et al. (2011 ) M,—L and M, — cr relations, which as noted 



above include BCG black hole masses, predict somewhat 
higher values of 1.1 x IQ^^Mq and 5 x lO^M©. It should 
be noted that A2261-BCG is more luminous and has a 
higher velocity dispersion than any galaxy that actually 
has a black hole mass measurement, thus the predictions 
are extrapolations. 
A different approach is to use the fundamental-plane 



of black hole activity ( Merloni et al.|2003 1 , which derives 
M, from a combination of 5 Ghz core radi o luminos- 
ity and 2—10 keV n uclear X-ray emission. Hlavacek 



Larrondo et al. (2012) use this methodology to estimate 
M. = 2.0+^;^ X iPW© for A2261-BCG. 

Unfortunately, estimates of black hole mass implied 
by the scale of the core, itself, in A2261-BGG will be 
even more uncertain extrapolations. Lauer et al. (200' 
presented a number of relationships between r^ and M, 
based on how the sample of core galaxies with real M, 
determinations was analyzed. If r^ and M, are fitted 



symmetrically, then Lauer et al. ( 2007a ) find 



which predicts M, ~ 4 x WMq for A2261-BCG. If on 
the oth er hand M, is trea ted as the independent variable. 



then Lauer et al. 



(2007a) find r^ cx M,^-^, which predicts 



a more modest M, ~ 7 x lO^M©. The large difference 
between the two estimates is due to the small number 



of systems contributing to the Lauer et al. (2007a) r^ — 
M, relationships and the large scatter between the two 
parameters. 

An alternate approach advocated in several papers 
is to relate estimates of the "mass-deficit," Md, that 
is the mass in stars that was ejected to generate the 
core, to M,. Unfortunately, this methodology at present 
also leads to large uncertainties in t he estimated M, for 
A2261-BCG. [Kormendy fc Bender] (2009) derive a re- 
lationship between Ld and M, for galaxies of the form 
Ld cx Af, , where Ld is the observed starlight rather than 
inferred mass deficit, and is measured by integrating the 
central difference between the observed surface bright- 
ness profile and an inward extrapolation of a Sersic-law 
fitted to the envelope of the galaxy. 

The best fit of a Sersic law to the envelope of A2261- 
BCG has n = 4.1, or essentially the classic r^/^-law form. 
An r^/^-law fitted to the envelope gives My » —20.8 
implying A/» ^ 1.1 x lO^^Af©, using the Kormendy & 



Bender ( 2009 ) relation. This later Ld is identical to the 
My Ki —20.8 derived by the sir nple estimate of "co re 
luminosity," = irl^r'^ 



Lauer et al. 



(2007a I as 



used by 

a proxy for Ld that avoids the need tor the caretul seiec- 
tion and evaluation of a "pre-scouring" reference surface 
brightness profile. 

4.2. Ejection of the Central Black Hole 

One intriguing possibility is that the large flat core 
has resulted from the central ejection of its nuclear black 
hole. In this scenario, the galaxy most likely had a large 
core to begin with, which would have been enlarged by 
the ejection of the central black hole. The large core 
would have been generated by a smaller black hole than 
that directly implied by the presumption that scouring 
b y the binary b lack hole did all the work. 



Merritt et al.| (2004) and .Boylan-Kolchin et al. ( |2004 1 
show that a core can be generated directly in a power- 
law" galaxy (a system that initially has a steep central 
cusp) when the components of a binary black hole merge 
and the remnant is ejecte d by the asymmetric emission 
of gravitational radiation. Gualandris & Merritt (20081 
simulated this scenario, finding that that ejection can 
substantially enlarge a pre-existing core, leading to the 
inference of an exceptionally large mass deficit. Given 
the high luminosity of A2261-BCG, it is extremely likely 
that it would have a large core prior to any ejection of 
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Figure 7. The relationship between cusp-radius and total galaxy luminosity (V-band Vega-based). The galaxy sample plotted was 
assembled in [Lauer ct al. (2007a) from a variet y of sources (the figu re is adopted from Figure 5 in that paper). The BCGs in particular 
come from the [Lame et al.| (|2U02J sample. The Lau er et al.| ||2007a[ r-y — L relationship (also given in equation H is plotted; the dotted 
lines indicate iitr scatter about the mean relationship. A2ZDi-hSUG is plotted at the top, clearly has a cusp-raoms larger than all other 
galaxies in the sample. The large core in the MS0735. 6-1- 7421 BCG discovered by |McNamara et al.| | j2009( l is also plotted for comparison. 

Real proof of the ejected black hole hypothesis would 
be to find direct evi dence of the ejected black hole, itself. 
Merritt et al. ( 2009 1 show that the ejected black hole car- 
ries along with it a "cloak" of stars that had previously 
been closely bound to it. The ejected black hole and asso- 
ciated stars would somewhat resemble a globular cluster 
or ultra-compact dwarf galaxy. The key diagnostic for 
such a system would be the extremely high velocity dis- 
persion of the tightly bound stars. An obvious question 
then for the present case is if any of the four sources 
proximal to the core might be such an object. 

An important point is that the stellar "cloak" is most 
likely to be considerably less massiv e than the black hole . 
For a black hole ejected from a core, Merritt et al. (2009 1 



a central black hole. Gualandris & Merritt (2008) also 
show that the central stellar density prohies can become 
extremely flat after ejection of the black hole, similar to 
what is seen in Figure [6] 

The physical displacement of the core from the enve- 
lope center, noted above, is indicative of a local dynam- 
ical disturbance. This is likely to be a relatively recent 
event. An ejected black hole would be trailed by a strong 
dynamical-friction wake as it leaves the core. In effect, 
it would "pluck" the core, which would then oscillate for 
a few crossing times, tc ^ r^/a or ~ 10^ yr for A2261- 
BCG. At the same time, this need not mean that the 
ejection of the black hole, itself, needs to be as recent. If 
the ejected black hole remains on a radial orbit that peri- 
odically returns to the core, then the disturbance might 
only be due to a recent passage (Merritt, private com- 
munication) . 
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find 

Mb 
M. 



5/2 



Vk 



10-^ km s' 



-5/2 



(5) 

where, Mb is the mass of stars bound to the black hole, 
a is the velocity dispersion of the galaxy, and Vk is 
the "kick" velocity with which the newly merged binary 
black hole is ejected from the center of the galaxy. In 
the case of A2261-BCG with a = 387 km s~^, the coeffi- 
cient in the above equation becomes 0.1. Whi le we have 
ech oed th e nominal Vk — 1000 km s^^ from |Merritt et 
al. ( 2009[ ), the real value depends on the unknown mass 
ratio and spins of the two black holes within the merging 
binary and can range from essentially zero to a few thou- 
sand km s"^. The only way to constrain Vk better would 
be to identify the remnant and measure its line-of-sight 
v elocity. 



Merritt et al. ( 2009 ) also show that the effective radius 
01 the cloak is also related to the kick velocity. For a 
black hole initially in a core, 

M, 



Re ~ 43 



lOiOM 







Vk 



103 km s' 



pc. 



(6) 



Given the expected scale and mass of the stellar cloak, 
the best candidate for an ejected black hole would be the 
least luminous of the four sources, the unresolved point 
source south of the core. The most luminous of the four 
sources has an implied luminosity of Ly 4 x 10^ Lq, 
and a substantially larger physical extent than the num- 
ber given above. When converted to a likely mass, 
the source carries a substantial fraction of the putative 
10^'^Mq black hole, thus not matching the expectation 
that Mb << M,. Morphologically, the source also resem- 
bles any number of other small galaxies visible within the 
neighborhood of the BCG. The two less luminous paired- 
sources, also north of the core, are more consistent with 
a possible cloak, but even their compact sizes are large 
for the expectations enumerated above. 

Of course if Vk is large, and the ejection happened long 
ago, then the remnant would be unlikely to be within the 
core, and the candidate list would have to be extended 
to sources at much larger distances from the center of 
the BCG. The time scale for exiting the core is simply 
r^/Vk, which is likely to be substantially shorter than 
tc- Conversely, a small kick favors larger and more ex- 
tended cloaks, and increases the likelihood that the rem- 
nant would remain close to the core in the event that 
the kick was insufficient to unbind the black hole from 
the BCG, itself. If the offset core is indeed due to an 
ejected black hole, then the short lifetime of the offset 
would imply that the remnant should be relatively close 

by. 

5. A2261-BCG AND THE FORMATION OF CORES 

There is a rich tradition in observational astrophysics 
of using the extreme member of an ensemble to under- 
stand the origin of the ensemble over all. The role of 
A2261-BCG as a witness to the mechanisms that form 
cores hinges on whether its core is "normal" or not. 
A2261-BCG bests its closest rival, NGC 6166, for having 
the largest core by a factor of two. More generally, the 
core of A2261-BCG stands out from all galaxies in the 



subset of 57 BCGs in the Lauer et al. (2007a) compos- 
ite sample of early type galaxies imaged wit h HST (the 



20021. Still, the 


the 

A Of" 


Lauer et al. 

m — ui '1 ' „ i„J 



3.2 kpc does not fall far above 
— L relationship, given that 



— even as compared to other BCGs. A2261-BCG might 
be expected to harbor a ~ IQ^^Mq black hole, based on 
its location within the M, — L and M, — a relations, thus 
the formation of core in this galaxy would a priori be ex- 
pected to a limiting case for the formation of the core by 
core scouring, the standard hypothesis. 

The core of A2261-BCG was not at first noted for its 
size, as much as its unusual appearance, however. Cores 
generally have at least a weak central surface brightness 
cusp, but in A2261-BCG there is no obvious center to the 
core. Its central stellar density profile is perfectly fiat, 
or perhaps even depressed at the center. While there 
is a radio source coincident with the core, there is no 
optical AGN counterpart or central nuclear star cluster 

— nothing suggests that A2261-BCG is hosting anything 
like a ~ lO^^M© black hole at its center. If binary black 
holes do scour out cores, then on occasion they must 
also merge and be ejected from the core, thus causing 
it to rebound and expand to an even larger size than it 
possessed at the completion of scouring. An attractive 
feature of this scenario is that it may account for large 
cores that may be difficult to exp lain by scouring alone. 
The simulations of Merritt et al. (2004) , 'Boylan-Kolchin| 
et al.| ( j2004J , and .Gualandris fc Merritt, ( ,2008) suggesF 
that the core of A2261-BCG matches the expectations of 
what a core that has ejected its central black hole looks 
like. 

There is a caveat, howe ver. Cores w i th eje cted black 
holes may be ephemeral. Faber et al. (1997) note that 
a central black hole can act as "guardian," "protecting" 
an existing core from being infilled by the central can- 
nibalism of less luminous, but centrally denser galax- 
ies. A2261-BCG lives in a rich environment. Figure 
[6] shows that the central mass density of the core is 
< 0.1 Mq pc-3^ which is extremely diffuse in comparison 
to the denser cores of less luminous galaxies. Without a 
central black hole in the core of A2261-BCG, the nuclei 
of galaxies merging wi th the BCG would readily settle 



intact into its center ( Holley-Bockelmann & Richstone 



2000 1 . The four sources projected against the core, are a 



reminder that the galaxy may cannibalize its neighbors, 
even if the sources in question, may themselves not be 
at risk for this. Again, if A2261-BCG did lose its central 
black hole, this may have been a relatively recent event. 

In the end, how we use A2261-BCG to test the the- 
ory of core formation hinges on whether or not it now 
hosts a supermassive black hole at its center. The radius 
interior to which a black hole dominates the stellar dy- 
namics is r, = GM,/a^ . For a ~ lO^^M© black hole and 
a = 387 km s~^, r, « 300 pc or w 0;'08 for A2261-BCG. 
This scale is readily accessible using adaptive optics in 
the near-IR on lOm-class telescopes, with the caveat that 
the extremely low surface brightness interior to the core 
will demand long exposures, even with such large aper- 
tures. The detection of a central M, ~ IO^^Mq would 
suggest that the core is simply an extreme example of the 
scouring mechanism. Conversely, a demonstration that 
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A2261-BCG lacks a suitably massive black hole, or the 
possible demonstration that one of the sources near the 
core is a cloaked black hole could establish that the large 
core indeed was generated by the ejection of its central 
black hole. An improved position for the central radio 
source might also determine if a black hole truly lies at 
the center of the core, or instead is associated with one 
of the knots. 
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